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Abstract 

We derive strong bounds on any possible large-scale spatial variation in 
the values of physical ’constants’ whose space-time evolution is driven by 
a scalar held. These limits are imposed by the isotropy of the microwave 
background on large angular scales in theories which describe space and 
time variations in the fine structure constant, a, the electron-proton mass 
ratio, n, and the Newtonian gravitational constant, G. Large-scale spatial 
fluctuations in the fine structure constant are bounded by 5a/a < 2 x 10~ 9 
and 5a/a < 1.2 x 10 -8 in the BSBM and VSL theories respectively, 
fluctuations in the electron-proton mass ratio by Sfi/fj, < 9 x 10 -5 in 
the BM theory and fluctuations in G by 8G/G < 3.6 x 10 _1 ° in Brans- 
Dicke theory. These derived bounds are significantly stronger than any 
obtainable by direct observations of astrophysical objects at the present 
time. 

PACS Nos: 98.80.Es, 98.80.Bp, 98.80.Cq 


1 Introduction 

The recent resurgence of interest in the possible slow variation of some tra¬ 
ditional ’constants’ of nature has focussed almost exclusively upon their time 
variation. This was led inspired to a considerable extent by the capability of new 
astronomical instruments to measure spectral lines created by the light from dis¬ 
tant quasar to very high precision. The quasar data analysed in refs. (Tj using 
the new many-multiplet method consists of three separate samples of Keck- 
Hires observations which combine to give a data set of 128 objects at redshifts 
0.5 < z < 3. The many-multiplet technique finds that their absorption spectra 
are consistent with a shift in the value of the fine structure constant between 
these redshifts and the present of Aa/a = [ a(z ) — a]/a = —0.57 ± 0.10 x 10 -5 , 
where a = a(0) is the present value of the fine structure constant. Exten¬ 
sive analysis has yet to find a selection effect that can explain the sense and 
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magnitude of the relativistic line-shifts underpinning these deductions. Fur¬ 
ther observational studies have been published in refs. -2] using a different but 
smaller data set of 23 absorption systems in front of 23 VLT-UVES quasars 
at 0.4 < z < 2.3 and have been analysed using an approximate form of the 
many-multiplet analysis techniques introduced in refs. JJ. They obtained 
A a/a = —0.06 ± 0.06 x 10 -5 ; a figure that disagrees with the results of refs. 
^]. However, reanalysis is needed in order to understand the accuracy being 
claimed. Other observational studies of lower sensitivity have also been made 
using OIII emission lines of galaxies and quasars. The analysis of data sets of 42 
and 165 quasars from the SDSS gave the constraints A a/a = 0.51 ± 1.26 x 10~ 4 
and Aa/a = 1.2 ± 0.7 x 10 -4 respectively for objects in the redshift range 
0.16 < 2 < 0.8 0. Observations of a single quasar absorption system at z = 1.15 
by Quast et al 0] gave Aa/a = —0.1 ±1.7 x 10 -6 , and observations of an absorp¬ 
tion system at z = 1.839 by Levshakov et al 0 gave Aa/a = 2.4±3.8 x 10 6 . A 
preliminary analysis of constraints derived from the study of the OH microwave 
transition from a quasar at z = 0.2467, a method proposed by Darling {B|, has 
given A a/a = 0.51 ± 1.26 x 10 -4 , j7j. A comparison of redshifts measured using 
molecules and atomic hydrogen in two cloud systems by Drinkwater et al |5j at 
2 = 0.25 and 2 = 0.68 gave a bound of Aa/a < 5 x 10~ 6 and an upper bound 
on spatial variations of 5a/a < 3 x 10 -6 over 3 Gpc at these redshifts; bounds 
on spatial variation of similar order arise from the results of ref. [1] because of 
the wide distribution of the target absorption systems on the sky. 

New observational studies sensitive to small variations in the electron-proton 
mass ratio, fi = m e /m p , at high redshift have also been reported 0Q331III], along 
with a new restriction on possible time variation of the Newtonian gravitation 
constant, G, in the solar system by the Cassini mission m A range of other 
astronomical and geophysical constraints which might limit possible changes in 
a have also been re-evaluated nanus]. 

Theories which can handle space-time variations in a and /r self-consistently 
(as opposed to simply ’writing in’ a time-varying constant into the equations 
in which it is truly a constant) have only recently been developed to comple¬ 
ment the Brans-Dicke gravity theory m which has been available to study the 
cosmological consequence of time-variations in G. The BSBM theory mm 
is a self-consistent extension of general relativity which incorporates space-time 
variations in the a, varying-speed-of-light (VSL) theories mi Emm m provide 
other ways of effecting variations in a and other gauge couplings, and the new 
theory of Barrow and Magueijo m allows space-time variations of /i to be stud¬ 
ied. All of these theories model the variation of a traditional constant by means 
of a scalar field which obeys a conservation equation derived from the variation 
of an action. The time variations of a, fi, and G that are permitted by these 
theories have been investigated in varying degrees of detail. They differ in one 
respect. The variations of the scalar fields carrying variations in a or /i do not 
have significant effects upon the expansion dynamics of the universe: the latter 
remains well described by the usual general relativistic FRW universe contain¬ 
ing the appropriate matter source. However, in the case of Brans-Dicke theory 
the changes in the associated scalar field do affect the expansion dynamics of 
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the universe and the FRW models are changed into new solutions that approach 
those of general relativity only in the limit that the space-time variation of G 
tends to zero. 

In this paper we will show how a simple treatment of these four scalar- 
tensor theories for varying a , p, and G allows us to predict and bound the 
magnitude of the spatial variations expected in these constants in the universe 
on extragalactic scales by using the observed temperature isotropy of microwave 
background radiation on large angular scales. We note that, in the past, bounds 
on spatial variation have been discussed by Tubbs and Wolfe [23, and Pagel, 
m who addressed this question at a time when the large-scale uniformity of 
the universe was a far greater mystery than it is in today’s post-inflationary era. 
These papers stressed that the values of combinations of physical constants that 
were found to be the same to high precision when deduced from the spectra of 
objects were so far apart on the sky that they could not have been in causal 
contact during the history of the universe prior to the emission of their light. 
Today, we expect a high degree of coherence within the whole of the visible 
universe because it may have evolved from the inflation of a single causally 
coherent domain. However, even if that were the case, if constants like a, p, 
and G are actually space-time variables and possess small quantum statistical 
fluctuations at the time of inflation then they may have a predictable (and even 
observable) spectrum of inhomogeneous variations today. A particular example 
is given by the chaotic inflationary universe in a Brans-Dicke theory of gravity 
EBBED, which gives rise to a spectrum of spatial fluctuations in the value of G 
as well as in the density of matter. 

In section 2 we shall describe four self-consistent theories of varying constants 
and in section 3 show how we can use the isotropy of the microwave background 
in conjunction with the predictions of these theories to derive bounds on the 
allowed spatial variations in these constants, before summarising our results in 
section 4. 


2 Four representative theories 

We will consider four representative scalar theories that are of particular interest 
given the current observational situation. It will be clear that these theories have 
an analogous structure. In each case the conservation of energy and momentum 
for the scalar field provides a wave equation of the form 

Dtp = \f((p)L(p,p), (1) 

where ip is a scalar field associated with the variation of some ’constant’ C via a 
relation C = f(p), A is a dimensionless measure of the strength of the space-time 
variation of C, f{p) is a function determined by the definition of p, and L(p,p) 
is some linear combination of the density, p, and pressure, p, of the matter that 
is coupled to the field p and f(p) — 1 for small p. At a given cosmic time, for 
small changes in p, this equation describes small spatial variations in C by a 
Poisson equation of the form 
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2.1 BSBM varying-a theory 

A simple theory with time varying a was first formulated by Bekenstein El as 
a generalisation of Maxwell’s equations but ignoring the consequences for the 
gravitational Held equations. Recently, this theory has been extended m to 
include the coupling to the gravitational sector and some of its general cosmo¬ 
logical consequences have been analysed. 

Variations in the fine structure ’constant’ are driven explicitly by variations 
in the electron charge, e, and the fine structure ’constant’ is given by 


where the the scalar ?/; field obeys an equation of motion of the form 0: 

Oip = e~ M p em , (3) 

UJi 

where ui is a dimensional constant which couples the kinetic energy of the i/j 
field to gravity and p em is the density of matter that carries electromagnetic 
charge. If we write £ = p em /p where p m is the total matter density then 

E 2 - B 2 
^ ~ E 2 + B 2 

and its sign (—1 < ( < 1) depends on whether the dominant form of (dark) 
matter is dominated by electrostatic ( E 2 ) or magnetic (B 2 ) energy (for further 
discussion see refs. mmmmm For the scalar field, we have the 
propagation equation, 

U* = - — e-^Cp m ( 4 ) 

CUi 

and for small variations in ijj and a this is well approximated by 

- —(Pm ( 5 ) 

U>1 

Some conclusions can be drawn from the study of the simple BSBM models 
with C < 0, [El- These models give a good fit to the varying a implied by the 
quasar data of refs. There is just a single parameter to fit to the data and 
this is given by the choice 


c_ 

Ull 


(2 ± 1 ) x 10" 4 


( 6 ) 


We shall use this as a conservative bound in what follows. Tighter observational 
limits will only strengthen our conclusions. 
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2.2 VSL theories 


In one ’covariant’ version of these theories ED ED ED variations in a (and all 
other gauge couplings a*) are driven by a scalar field y that drives explicit vari¬ 
ations in the speed of light and couples to all the matter fields in the lagrangian, 
not just the electromagnetically-coupled matter, with cc* = exp[Qy], where Q is 
a numerical constant. The structure of covariant VSL is analogous to the BSBM 
theory and for small variations in the fine structure constant, exp[\] ~ 1 , we 
have 


Ox-~—Pm (7) 

U) 2 

where u >2 is a coupling constant. The observational data of Webb et al [Tj are 
fitted by 


Q_ 

UJ2 


= 8±4 x 1CT 4 


which we use as the observational bound on the coupling. There is no variation 
of a in the limit that ^ —*> 0. Another edition m of a VSL theory has variation 
only in the electromagnetic coupling, a, and variations are driven only by the 
pressure of matter: 


□x - -47 rGujp(p) (8) 

for some new coupling constant to. We will just examine the covariant VSL 
theory, eq. 0 in what follows. 

2.3 BM varying-/i theory 

The theory recently devised by Barrow and Magueijo m describes a varying 
electron-proton mass ratio, /r, via a changing electron mass which is driven by 
a scalar field, <f >, defined by 


m e = moe^ 


where <f> obeys 

D ^_ mo(n e -n p ) ^ 
w 3 

Here, n e and n p are the electron and proton number densities, and 0*3 is a 
dimensional coupling constant. In the case of small variations (e^ ~ 1) this is 
well approximated by 


nii>~nn~-£z (9) 

UJ3 

and observational bounds on the time-variation of fj, at high redshift u impose 
a weak bound of 
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Guj 3 > 0.2 


The 0 J 3 —> 00 limit is that of constant p. 

2.4 Brans-Dicke gravity theory 

The Brans-Dicke theory m generalises Einstein’s general theory of relativity 
to incorporate a space-time variation in the Newtonian ’constant’ G by means 
of a Brans-Dicke scalar field $ oc G _1 which obeys a conservation equation of 
the form 


d$= Mp-3 P ) i 

3 + 2 u>bd 


( 10 ) 


where p and p denote the total density and pressure of matter respectively, u>bd is 
the dimensionless Brans-Dicke parameter and general relativity (G = constant) 
is approached as u>bd —> 00 . The current observational lower bounds on the 
allowed time variation of G give udbd > 500 from a variety of local gravitational 
tests (see JT3] for a review). But the strongest constraint to date is derived from 
observations of the time delay of signals from the Cassini spacecraft as it passes 
behind the Sun. These considerations led Bertotti, less and Tortora m, after 
a complicated data analysis process, to claim that 


uj bd > 40000 (2a). 


( 11 ) 


This theory differs from the three listed above in that small variations in $ 
have a significant effect upon the expansion dynamics of the universe because 
these variations control the strength of gravity. These variations can be seen 
explicitly by writing down the Friedmann equation for the expansion scale factor 
a(t) in the case of zero spatial curvature: 


a?_ _ 87t <j> u>bd 4> 2 

a 2 ~ 34> P <f> + 6 $ 2 ' 

For the case of dust (p = 0) there are simple power-law solutions: 


( 12 ) 


= £[2+2tc>b c 2]/[4+3 Ubd] 


(13) 


$(t) = $ 0 t 2 /[ 4 +3^w] ( 14 ) 

and we note that we recover the usual Einstein-de Sitter cosmology with a = t 2, :i 
as u>bd —> 00 . These solutions solve eq. m exactly in the p = 0 case m 


□<i> 


8 np 

3 + 2 u>bd 


(15) 
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3 Large-scale inhomogeneity in a, /i, and G 

The four theories that we have introduced have a characteristic structure in 
which changes in cr, /z, and G are driven by different parts of the density of 
matter in the universe. This means that inhonrogeneity in the matter content 
of the universe is coupled to inhomogeneities in the ’constants’ a, /z, and G. 
If we ignore the non-Machian mode that is not driven by the matter fields 
(ie the complementary function arising from the solution of D(scalar) = 0) 
because it falls off rapidly in time and becomes negligible at late times in the 
universe, then we can estimate the allowed large-scale spatial inhomogeneity 
of a,/z, and G in terms of observable quantities. The inhomogeneity in the 
’constants’ requires inhomogeneity in the driving matter perturbations and their 
associated gravitational potential fluctuations. The magnitude of the latter 
is observationally constrained by the temperature isotropy of the microwave 
background on large angular scales and leads to an upper bound on the possible 
inhomogeneity in the values of ’constants’. We shall ignore the acceleration of 
the universe which began recently at z ~ 0.3. Its inclusion leads to a very small 
change in the final results, no larger than the uncertainty in other parameters. 

3.1 BSBM inhomogeneity in a 

On a hypersurface of constant comoving proper time, t , eq. © gives 


V 2 W) ^ ^ 5 P] 


So, noting that QirGpt 2 — 1, any inhonrogeneity in a over a scale L is linked 
to inhomogeneity in the density of electromagnetically-coupled matter, Sp em , 
and/or inhomogeneity in the electromagnetic quality of the dark matter, 5(^ ern , 


by 1 



However, we expect that the inlromogeneity in the electromagnetically coupled 
matter can be written in terms of the inhonrogeneity of the total matter density 
on scale L in terms of some biasing parameter (3 which will not depart too 
greatly from being 0(1), so we put 


^ Pern ^ p ^Pm. 


era 


Pem Pm 


Hence, we have 


1 Note that V 2 i/> —► (V'M — ' l P( x )) where is the average value of ^ in a spherical 


region of radius L as L —> 0 and is defined by i^{x) = + r)d?rj. 
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5a £ 

— ~ 0.3— 

a lo i 


m 

c 



but we note that the gravitational potential perturbations (5&n/&n) in the 
space-time metric created by 5p m are given by Poisson’s equation (V 2 <f>jv = 
4 tt Gp m ) as 


(-Y 

P m \t ) 

and these fluctuations produce temperature fluctuations in the microwave back¬ 
ground radiation, A T/T ~ 5^n/^n — 2 x 10~ 5 , on the corresponding angular 
scales, m so that 


5<$>]\[ 


0.3 


Hv 



~6x 10“ 5 


on large angular scales (9 > 10°), [33]. We will assume, in accord with obser¬ 
vations, that the gravitational potential fluctuations are scale-independent to a 
good approximation 2 . 


Hence, if we take the best fit to 


_C_ 

Ui 


2 x 10 4 from eq. © the observations 


of ref. [2 and we assume that the spatial variations in the electromagnetic 
composition of the matter in the universe are approximately proportional to 
the variations in the matter density, with 


HO ^ Hem 

~ ? ?-> 

Pern 


where p ~ 0(1), then the spatial fluctuations in a and the required microwave 
background temperature anisotropies in these theories are simply related by 

— ~ 0.9—/3(1 +i,)^-2x 10 -9 /3(l + rj). (16) 

a u) i T 


Hence, for the two cases of small and large spatial variations in £, respec¬ 
tively, we have 


5a 

2 x 10 

a 


5a 

2 x 10 

a 



~ 9 l 

—9 r. 


2 Note that the smallness of <5 $jv/^jv is the justification for the so called ’cosmological 
principle’ and the use of the Friedmann metric. The smallness of 5p rn /p rn is unnecessary 1)51 . 













These bounds on large-scale spatial variations of the fine structure ’constant’ 
are extremely strong and we have obtained them by assuming there is a bound 
on the level of time-variation consistent with the observations of ref. ,I|. For 
comparison, the sensitivity of direct searches for variations in a which compare 
observations of different quasar absorption spectra around the sky is only about 
da/a < 0(1CT 6 ). Recall that, in contrast, the time variation of the fine struc¬ 
ture constant is far more strongly constrained by the quasar absorption-system 
data than by the microwave background effects on small scales m, 

3.1.1 VSL inhomogeneity in a 

A similar argument to that used for the BSBM case can be applied to the 
covariant VSL theory and leads to the result that the allowed spatial variation in 
a is again bounded by the temperature anisotropy in the microwave background 

by 


5a 

a 


0 . 3 - 2 . (TV 
^2 Pm \t J 


0.9 Q- — < 1.2 X 1CT 8 . 
U)2 T 


A similar bound could be deduced for the allowed spatial variations in all the 
gauge couplings, Sat/ati in this theory. Again, this bound is far stronger than 
can be achieved by direct spectroscopic studies of quasars and other astrophys- 
ical systems at low redshift. 


3.2 BM inhomogeneity in /i 

If we apply this argument to possible spatial variations in the electron-proton 
mass ratio, fi 1 described above then a similar argument leads to an expression 
for the inhomogeneity in the electron-proton mass ratio, /i, of the form 


6fi 

P 


0.3 I" Sp e (L_V 
Guj 3 p e \t) 


0 . 9/3 AT 
Goj 3 T 


< 9 x 10“ 5 (3 


where we have assumed that the inhomogeneity in the electron density is ap¬ 
proximately proportional to that in the matter distribution: 

$Pe = p 5 Pm 
Pe Pm 

In this case the numerical bounds on the allowed variation are much weaker than 
for 5a/a. This is a reflection of the weak constraints that exist on time-varying 
/i in these theories m because of the different time-evolution of a and /r during 
the radiation era. 


3.3 Brans-Dicke inhomogeneity in G 
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The analysis of the level of inhomogeneity expected in a Brans-Dicke universe is 
slightly different because time variations in <f> ~ G~ l determine the expansion 
dynamics. For the ’Machian’ solutions m-m that are the attractors at late 
times we have 87 jp ~ and so inhomogeneity in the Brans-Dicke held (<5<f> ^ 
0) is linked to inhomogeneity in the gravitation ’constant’ (SG ^ 0) and in the 
total matter density (Sp ^ 0 ) by 

SG _ 5$ ~ 10.8 S Pm (L\ 2 {u bd + l ) 2 10.8 (co bd +1 ) 2 AT 

G 3 + 2ubd Pm V t ) (4 + 3ubd ) 2 3 + 2 ivbd (4 + 3wm ) 2 T 

Hence, for large values of wm, as observations of the time-variation of G require, 
this simplifies to 


SG 

~G 


- x 10 " 5 < 3.6 x 10 " 10 

5 tobd 


if we use the Cassini bounds on u>bd ■ If we replace the Cassini bound by an obser¬ 
vational bound of u>bd > 1000 from other solar-system constraints this weakens 
the bound by a factor of 40. In either case the bound on spatial variations is 
extremely strong. It results from a combination of the microwave background 
anisotropy limits on the density perturbations which drive variations in G and 
the intrinsic weakness of the $ coupling. Notice that, just as fluctuations in 
the microwave background temperature are known far more accurately than 
the mean temperature itself, so the spatial fluctuations in G are limited to far 
greater accuracy than the value of G is known experimentally (see for example 

E3). 

We summarise the principal results obtained in Table 1. 


Theory 

Scalar field 

Scalar coupling 

Scalar equation 

Bound on inhomogeneity 

BSBM 

ip : a = e 2 ^ 


X 

Oil 

= 2+ 4 xl0 -4 

ap = 2p - e ~" 

Ai~X / 3 ?7 AT< 2 x 10" 9 

VSL 

X : a t = e Qx 


Q_ 

CO 2 

1 

O 

X 

H“ 1 

OO 

II 

^ CO 2 

12 x 10~ 8 

a C 02 T ~ 

BM 

<p : p oc 

Gw 3 > 0.2 

n<p= 

co 3 

9 x 10~ 5 

u Geo 3 T ~ 

BD 

$ : G oc 4) -1 

u) bd > 40000 

_ 3+2 ui hd _ 

, 3 A T T <3.6x 10" iU 

_ G _5 uJhd 1 ~ _ 


Table 1: Summary of the inhomogeneity levels allowed by the microwave- 
background temperature isotropy in the theories discussed in this paper. In 
column 2 the link between the scalar field and the ’constant’ is defined; in 
column 3 observational bounds on the dimensionless scalar coupling constant 
appearing in the theory and discussed in the text are given; the propagation 
equation for the scalar field is given in column 4 and in column 5 the bounds on 
any spatial variation in the associated constant imposed by the bound on the 
coupling and the isotropy of the microwave background (A T/T < 2 x 10 -5 ) are 
summarised. 


10 
























4 Discussion 


We have shown that the characteristic structure of scalar theories for the space- 
time variation of supposed constants of Nature enables us to use the observed 
bounds on the couplings in the theories obtained from observational searches 
for time variations in the associated ’constants’ and the microwave background 
isotropy to obtain very strong bounds on any spatial fluctuations in these con¬ 
stants over large astronomical scales. Since the scalar fields which carry the 
space-time variations of constants are driven by all (or part) of the matter con¬ 
tent of the universe there will always be density inhomogeneities present which 
drive inhomogeneities in the constants. However, the gravitational potential 
fluctuations associated with these large-scale density perturbations show up 
in the microwave background temperature anisotropy on large angular scales. 
Their amplitudes are therefore bounded above and this bound in combination 
with limits on the strength of the coupling of the scalar field leads to a series of 
very strong bounds on possible spatial fluctuations in the associated constants. 
We have calculated these bounds for four representative theories which are of 
current interest. The basic argument is of wider application to other dilaton 
theories and with small modifications the bounds can be modified to include 
the small changes that arise because of a non-constant curvature spectrum of 
density perturbations. They will be strengthened if the underlying couplings in 
these theories are strengthened by future or present observational studies. 

Detailed consideration of the behaviour of fluctuations on smaller scales 
would lead to a different range of constraints. We note that the limits pre¬ 
sented here apply to large scales and do not apply to the possible variations 
in the values of ’constants’ that can arise because of the non-linear evolution 
of cosmic overdensities into clusters, galaxies and planetary systems. These 
are not bounded by the isotropy of the microwave background and still permit 
significant spatial variations to arise on small scales in the local universe M- 
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